Observation of coupled normal heat and matter transport in a simple model system 
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We introduce the first simple mechanical system that 
shows fully realistic transport behavior while still being ex- 
actly solvable at the level of equilibrium statistical mechanics. 
The system under consideration is a Lorentz gas with fixed 
freely-rotating circular scatterers which scatter point particles 
via perfectly rough collisions. Upon imposing either a tem- 
perature gradient and/or a chemical potential gradient, a sta- 
tionary state is attained for which local thermal equilibrium 
holds. Transport in this system is normal, in the sense that 
the transport coefficients which characterize the flow of heat 
and matter are finite in the thermodynamic limit. Moreover, 
the two flows are non-trivially coupled, satisfying Onsager's 
reciprocity relations to within numerical accuracy. 
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The phenomenological equations of irreversible ther- 
modynamics, which describes diffusion, heat conduction, 
viscosity among a host of other phenomena, are of great 
significance in many fields. To link these phenomena to 
the underlying microscopic dynamics is, however, highly 
non trivial and the results so far have not been conclu- 
sive. From a mathematically rigorous point of view, very 
few results have been obtained From the numerical 
point of view, the systems considered thus far are either 
too complicated to shed much light upon the problem, 
or fail to reproduce the macroscopic phenomenology. In- 
deed, to our knowledge, the validity of Fourier's law has 
been proven analytically only for a specific model in the 
limit of infinite dilution with finite mean free path 0,0. 
There have also been attempts to link transport phenom- 
ena to the chaotic properties of the underlying classical 
dynamics [0, and a connection between the rate of en- 
tropy production and the rate of contraction of phase 
space volume has also been pointed out ||. Further re- 
sults have also been obtained concerning energy trans- 
port in the Lorentz gas [||. However, as was shown in [Q, 
this last system does not satisfy local thermal equilibrium 
(LTE), which is a fundamental assumption of irreversible 
thermodynamics. It is therefore not clear what precise 
meaning should be attached to Fourier's law in this situa- 
tion. Numerically, attempts have been made, on one side, 
through the simulation of realistic many-body systems 
satisfying a thermostatted (non Hamiltonian) dynamics 
||. These simulations have indeed been able to repro- 
duce non-trivial transport phenomena. However, such 
studies do not provide a detailed understanding of the 
microscopic processes involved due to the excessive com- 



plexity of the system. The other numerical approach in- 
volves the study of transport in "simple systems" . Exam- 
ples of these include: chains of anharmonic oscillators 
for which Fourier's law does not hold; Lorentz gases [p|, 
which, as mentioned above, do not satisfy LTE; and the 
so called ding-a-ling and ding-dong models Jl0|,[ll| , which 
indeed yield normal heat transport, but have geometric 
constraints that make even their equilibrium properties 
a complicated affair and cannot be generalized, retaining 
their simplicity, to obtain the coupled heat and matter 
transport commonly found in real systems. Under these 
circumstances, the minimal ingredients required in the 
microscopic physics of a system to attain normal heat 
conduction in low dimensions arc presently still under 
discussion Qj^fLf]. 

In this work we introduce the first simple reversible 
Hamiltonian model system whose steady state is well de- 
scribed by the hypothesis of LTE and supports both nor- 
mal heat and matter transport. These are non-trivially 
coupled, satisfying Onsager's reciprocity relations. Fur- 
thermore, this system has the advantage that its equilib- 
rium properties are trivial since its description reduces to 
that of ideal gases. This makes our model an ideal testing 
ground for theories linking microscopic mechanical prop- 
erties to transport phenomena. Moreover, this model can 
be easily modified to study other physical problems from 
a microscopic approach, such as transport in heteroge- 
neous systems |14 and the effects of applied fields. 

Our model consists of non-interacting point particles 
of mass m scattered by an array of freely-rotating circular 
scatterers of radius R with a finite moment of inertia O. 
The only relevant adimensional parameter characterizing 
the system is 
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which determines the amount of energy transfer between 
disks and particles in a collision. The scattering proceeds 
according to the rules characterizing perfectly rough col- 
lisions, designed to conserve energy and angular momen- 
tum. They are given by the following prescription: the 
normal velocity v n of the particle changes sign, whereas 
the tangential velocity vt of the particle and the angular 
velocity u of the disk are transformed as follows: 
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Rjj' = Ruj^ (v t - Rlo) . 

1 + 77 

These rules define a deterministic, time-reversible, and 
canonical transformation at each collision. For finite val- 
ues of 77, particles in the system may exchange energy 
among each other through the disks, even though they 
do not interact directly. It is precisely this disk-mediated 
energy exchange which permits the systems to reach ther- 
modynamical equilibrium without the necessity of any 
additional thermostat (a similar idea was proposed in 
|jl6f as a model of a deterministic thermostat). It should 
be noted that both the limits 77 — > and 77 — > 00 are 
exceptional: In the former case the particle dynamics re- 
duces to that of normal elastic scattering by a hard disk, 
which is the usual Lorentz gas model; in the latter, the 
state of rotation of the discs is unaffected by the colli- 
sions. Thus, in both limits the energy-mediating effect 
of the disks is suppressed and thermodynamical equilib- 
rium is not reached. In the following, unless the contrary 
is explicitly stated, we shall always be dealing with the 
case 77 = 1, since that is the value of 77 for which equili- 
bration is most efficient. We shall also always set m and 
R to one. 

A peculiarity of this system is that it is homogeneous, 
i.e. it has no proper energy scale, or, equivalently, no 
proper time scale. Thus all energies and temperatures 
reported in this work are rescaled to the lower nominal 
temperature of the baths. 

The geometric disposition of the scatterers is indicated 
in Fig. [j]. The centers of the scatterers are fixed on a 
triangular lattice, along a narrow channel with periodic 
boundary conditions in the vertical direction. At the 
ends, the walls are used to fix both the temperature T 
and the quantity fi/T, where /j is the chemical potential. 
This is achieved by ensuring that each particle that hits 
the right or the left wall is absorbed with probability 
1 — exp(—l/\v n \), where v n is the normal component to 
the wall of the particle velocity. Otherwise, the particle 
is reflected with a velocity chosen from 

P„K) = ^K|exp^-|^ , 

where v n and vt are, respectively, the normal and tan- 
gential components of the velocity with respect to the 
wall, and T is the nominal temperature of the wall. Ad- 
ditionally, the walls emit particles with velocities dis- 
tributed as in (j^), at a velocity dependent rate given 
by 7(1 — exp(— l/|?j ra |). Here, 7 determines the nomi- 
nal value of the chemical potential at the walls through 
fx = 7Tn(7/T 3 / 2 ). Details will be given in a forthcom- 
ing publication []15[ . At the level of equilibrium statisti- 
cal mechanics the system is an ideal gas mixture, with 
the disks having one degree of freedom and the particles 




FIG. 1. Schematic illustration of the scatterer geometry: 
the scatterers are disposed on a triangular array with finite 
horizon to avoid infinitely long trajectories. For matters of 
convenience, in this work the separation between scatterers 
has been set to have the critical horizon. Periodic boundary 
conditions are used in the vertical direction. To avoid spurious 
effects arising from multiple scattering off the same disk, we 
have put two disks on each vertical. To study the dependence 
with system size, the length L of the sample is varied. The 
quoted length is the number of disks. 

two. We have verified that our model reaches an equilib- 
rium state fulfilling the equipartition theorem in the mi- 
crocanonical, canonical and grandcanonical ensembles of 
equilibrium statistical mechanics. Further, the tempera- 
ture and chemical potential at which the system equili- 
brates coincides with the nominal values at the walls in 
the canonical and grand canonical cases. Finally, when 
the system is subjected to a weak temperature or chemi- 
cal potential gradient, it reaches a well-defined stationary 
state after a relaxation time that depends on the length 
of the system and on the gradients themselves. To show 
that LTE is indeed achieved in the stationary state at- 
tained by our system, we display in Fig. ^| a typical re- 
alization of the energy distribution function P x (e) of the 
particles as they cross a narrow slab centered around po- 
sition x for three different positions. At each position, 
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FIG. 2. Semilogarithmic plot of the particle's energy distri- 
bution P x (s) at different positions along a channel of length 
30 with an end-to-end temperature difference AT = 1/7 and 
A(fi/T) — —0.2002. The different curves correspond to a fit 
to a Boltzmann distribution for each position. From these fits, 
we obtain the temperatures: T = 1.0193 at x — 3 (circles), 
T = 1.0693 at x = 13 (squares) and T = 1.1335 at x = 27 
(triangles). The curves have been rescaled for clarity. 
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FIG. 3. Temperature profile for the simulation described in 
Fig. |^. The solid line corresponds to the particle temperature 
along the channel while the open circles are the temperature 
of the disks, averaged over 500 realizations. In the inset the 
particle temperature profile (solid line) is compared with the 
temperatures obtained from the fits to the Boltzmann distri- 
bution of -Pc(e) in Fig. H, (open squares). 



P x (e) is consistent with the Boltzmann distribution. Fur- 
ther, if we determine the temperature T{x) by a fit of 
the P x (s) to the Boltzmann distribution, the agreement 
with the expected linear temperature profile in the sys- 
tem is satisfactory and coincides with the average energy 
of the particles measured at that position, see the inset 
of Fig. ||. Thus, the identification of the average parti- 
cle energy with the local temperature is justified in this 
system. In Fig. || we show the temperature profiles of 
both particles and disks, which are found to be linear. 
The agreement between both values indicates that disks 
equilibrate with the particles locally along the channel. 
In the inset, the temperatures obtained in Fig. || for three 
different positions are compared with T(x), measured as 
the average energy. The agreement reinforces the con- 
clusion that the system establishes a LTE in its steady 
state. Note again that T{x) reaches the nominal tem- 
peratures specified by the walls at both ends. This is in 
contrast to what is observed in several other models for 
heat conduction previously proposed JTf|,p|JT^] . Moreover, 
in simulations in which we impose a temperature gradi- 
ent but the same value of \xjT on both walls, we find 
that p(x)/T(x) is constant along the channel, confirm- 
ing that the particles in the system can be described as a 
two-dimensional ideal gas which is locally at equilibrium. 

From the general theory of irreversible processes, to 
linear order the heat and particle currents J u and J p can 
be written as follows (see e.g. pTy): 

Ju = L UU X7— — L up V — , 
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FIG. 4. Size dependence of the heat and matter currents for 
simulations with a fixed temperature difference, AT = 1/7, 
and fi/T constant. The dotted lines corresponds to 1/L scal- 
ing. From the agreement it follows that the transport coeffi- 
cients do not depend on the size of the system. 



1 fJ. 
Jp = L pu V — — LppV — 



(4) 



and the Onsager reciprocity relations read in this case 
Lup = Lp U . (5) 

The central feature of the model is that its transport 
properties are normal, meaning that the various trans- 
port coefficients appearing in (^) do not depend on the 
length of system, and are, thus, well defined in the ther- 
modynamical limit. In Fig. 4 we show the dependence of 
the currents on the length L of the system, for a typical 
realization. The 1/L dependence observed confirms that 
transport is normal. 

In order to obtain the value of the coefficients in (Q), it 
is enough to perform two simulations: Fixing the value 
of VT and setting V(/z/T) = yields L uu and L pu while 
setting VT = and fixing V(/i/T) gives L up and L pp . 

We have performed simulations with temperature and 
chemical potential differences up to 20% of the minimal 
nominal values at the walls, and in all cases we have found 
normal transport consistent with (||). For example, in a 
simulation in which a temperature difference of AT = 
1/14 and L — 30 at fixed fi/T = —5.01, averaged over 
500 realizations, we found that after the steady state has 
been reached, both a heat current and a particle current 
were driven by the temperature gradient. The measured 
value for L pu in this case was 0.0054 ± 0.00013. This 
implies, through that a heat current will be driven 
by a gradient in /i/T at fixed T, as is indeed observed. 
In the complementary simulation with T constant, where 
Af = -0.06, we found L up = 0.0055 ± 0.00013, thus 
confirming (|^) to within our numerical accuracy. 

As a consistency check, we have also studied a "canon- 
ical" situation, in which we supressed absorption and 
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emission of particles at the walls. In this situation there 
is no flow of matter in the steady state. The relationship 
between heat flow and temperature gradient becomes 
J u = kVT, with k given by the following expression: 



T 2 L 



(6) 
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This relationship was found to hold to good accuracy, 
thus confirming the validity of our "grand canonical" sim- 
ulations by which the L's were evaluated. 

We argue that the coupling between the two currents is 
non-trivial in the following sense: in a "canonical" simu- 
lation, the simplest assumption for the dependence of the 
density on the position is that the orbit of each particle 
covers the sample uniformly. Then the local temperature 
merely determines the speed at which the orbit is being 
traversed. This would imply that in such a simulation, 
the particle density scales inversely with the average ve- 
locity, that is 



p{x)T(x) x / 2 = const . 



(J) 



In terms of the transport coefficients defined in (|4|), 
equation (^) is equivalent to L pu = (d+l)TL pp /2, where 
d is the dimension of the system. This relation corre- 
sponds to a system for which all transport arises from 
uncorrelated Markovian motion of the particles, as in the 
Knudsen gas yjj . 

However, equation (Q) does not hold in our system as 
we always find a systematic spatial variation in this quan- 
tity, the size and sign of which depend on the value of 
n. Thus, our system cannot be accurately described in 
this simple manner. The correct description is, as in all 
systems showing realistic transport properties, an open 
problem. 

In summary, we have introduced a reversible Hamilto- 
nian system to study transport, with the following de- 
sirable properties: Its statistical mechanics is trivial, re- 
ducing to that of ideal gases. In the steady state it is 
well described by the hypothesis of local thermal equilib- 
rium. Its transport properties are normal and it naturally 
supports coupled matter and heat transport for which 
Onsager's reciprocity relations are satisfied to within nu- 
merical accuracy. This coupling is non-trivial in that it 
does not reduce to a simple velocity scaling. In view 
of these properties, we believe that this model provides 
an ideal framework from which theories for the emer- 
gence of macroscopic transport in out-of-equilibrium con- 
ditions can be constructed and/or tested. Furthermore, 
the model opens the possibility to study other phenom- 
ena such as thermal junctions (and other kinds of inho- 
mogeneities) by changing the size or moment of inertia of 
the disks, Joule heating by applying an electric field, bro- 
ken time reversal symmetry by applying magnetic fields, 
etc jl8) . Research into some of these extensions is cur- 
rently underway [|15| . 
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